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•  The  performance  of  SOFC  with  non- 
uniform  potential  operation  was 
evaluated. 

•  Maximum  stack  temperature 
gradient  and  stack  temperature  were 
considered. 

•  Non-uniform  potential  SOFC  stacks 
with  a  cooler  show  higher  perfor¬ 
mance. 
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The  performance  and  the  operation  viability  of  solid  oxide  fuel  cells  (SOFC)  fueled  by  reformed  methane 
using  the  concept  of  non-uniform  potential  operation  were  evaluated  by  one-dimensional  finite  differ¬ 
ence  analysis.  Two  thermal  constraints  including  maximum  acceptable  stack  temperature  gradient  and 
maximum  stack  temperature  were  the  criteria  for  considering  the  operation  viability.  The  average  power 
density  was  used  as  a  performance  indicator.  In  addition,  two  major  process  parameters:  fuel  and  air 
inlet  temperatures  and  the  operating  voltage  were  examined  to  determine  the  suitable  values  for 
different  viabilities.  The  results  revealed  that  non-uniform  potential  SOFC  with  two  stack  sections 
operating  at  different  voltages  could  not  offer  a  significantly  higher  average  power  density  compared 
with  the  uniform  potential  SOFC  owing  to  thermal  constraint  limitations.  The  installation  of  coolers 
between  the  stacks  is  important  and  it  could  offer  12.8%  improvement  in  average  power  density. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Due  to  the  continuous  increase  in  fossil  fuel  price,  alternative 
developments  of  high-performance  and  environmental-friendly 
technologies  for  electricity  generation  are  desirable.  Fuel  cell 
technologies  have  been  of  interest  over  the  last  decade  owing  to 
higher  electrical  efficiency  compared  with  conventional  combus¬ 
tion  heat  engines  [1].  This  is  due  to  the  direct  conversion  of 
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chemical  energy  to  electrical  energy  via  electrochemical  reactions 
whilst  a  prior  conversion  of  the  chemical  energy  to  thermal  and 
mechanical  energies  takes  place  in  the  combustion  heat  engines 
along  with  greater  irreversible  losses.  Among  various  types  of  fuel 
cells,  the  solid  oxide  fuel  cell  (SOFC)  is  a  promising  technology  as  it 
offers  number  of  advantages:  (1)  high  operating  temperature 
which  enhances  the  rate  of  electrochemical  reactions  and  inhibits 
the  overpotential  losses,  (2)  improved  electrical  efficiency  as  the 
high-temperature  exhaust  gas  (1073-1273  I<)  is  utilized  to  produce 
additional  electricity  via  gas  turbine/engine  cycle  [1,2],  and  (3)  fuel 
sources  such  as  methane,  ethanol,  natural  gas  or  oil  derivatives  can 
be  used  without  pre-processing  constraints. 
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The  improvement  of  SOFC  performance  has  been  performed  using 
three  main  methods:  the  development  of  a  new  component  material 
to  offer  a  low  overpotential  loss  [3,4],  the  combination  of  SOFC  stack 
with  heat  engines  to  generate  additional  electricity  [5,6]  and  the 
optimization  of  operating  conditions  [7].  Among  the  optimization 
techniques,  concept  of  using  the  so-called  “non-uniform  potential 
operation”  is  of  interest  [8-11  ].  In  this  concept,  a  SOFC  stack  is  divided 
into  several  sections  and  operated  at  different  operating  voltages. 
Near  the  entrance  of  the  fuel  cell  stack,  where  H2  concentration  in  fuel 
gas  and  rate  of  electrochemical  reaction  are  high,  a  high  value  of  cell 
voltage  is  used  to  maximize  the  power  density.  Inversely,  near  the  exit 
of  the  stack,  low  cell  voltage  is  set  to  maintain  SOFC  power  density  at  a 
high  value  although  the  electrochemical  reaction  rate  is  low  due  to 
low  H2  concentration  of  the  fuel  gas.  Selimovic  and  Plasson  [8] 
investigated  the  non-uniform  potential  operation  in  networked- 
SOFC  stacks  combined  with  a  gas  turbine  cycle.  They  concluded  that 
an  increase  in  the  electrical  efficiency  of  5%  could  be  achieved  by 
improving  the  thermal  management.  Also,  Au  et  al.  [9]  found  that  0.6% 
improvement  in  the  electrical  efficiency  could  be  achieved  for  molten 
carbonate  fuel  cell  (MCFC)  when  the  concept  of  non-uniform  potential 
operation  was  applied  [9].  The  implementation  of  this  concept  in 
polymer  electrolyte  membrane  fuel  cell  (PEMFC)  was  also  examined 
[10].  The  results  indicated  that  PEMFC  implementing  the  concept  of 
non-uniform  potential  operation  offered  6.5%  improvement  in  the 
power  density  compared  with  the  conventional  PEMFC.  Vivanpatar- 
akij  et  al.  [11]  investigated  the  implementation  of  this  concept  by 
considering  power  density  as  a  performance  indicator  and  achieved 
9.2%  enhancement  in  the  power  density  by  implementing  a  non- 
uniform  potential  operation.  Although,  non-uniform  potential  oper¬ 
ation  improved  the  performance  of  SOFC,  further  operating  parame¬ 
ters  such  as  the  thermal  behavior  of  solid  part  inside  SOFC  stack 
implementing  this  concept  would  promote  even  more  the  perfor¬ 
mance  as  a  large  amount  of  heat  is  generated  due  to  the  irreversibility 
of  the  electrochemical  mechanism  at  higher  level  of  the  electrical 
generation.  Recently,  Guan  et  al.  [12]  developed  an  experimental 
method  to  measure  temperatures  inside  stacks  using  thin  K-type 
thermocouples  and  self-developed  CAS-I  sealing  materials.  The  re¬ 
sults  have  shown  that  temperature  distribution  of  unit  cells  was  non- 
uniform:  the  maximum  temperature  difference  between  inlet  and 
outlet  was  about  200  K.  The  lack  of  knowledge  in  the  selection  of  the 
operating  conditions  would  cause  high  temperature  gradients  along 
the  SOFC  stack  and  some  solid  parts  in  the  stack  would  be  damaged. 

In  this  study,  the  implementation  of  “non-uniform  potential 
operation”  concept  in  SOFC  fueled  by  pre-reformed  methane  was 
examined  using  one-dimensional  analysis  (1-D  analysis)  for  the 
investigation  of  stack  performance  and  thermal  behavior  along  the 
SOFC  stack.  The  calculation  using  1-D  analysis  could  well  represent 
the  actual  SOFC  stack  behavior  as  reported  by  Sorrentino  et  al.  [13]. 
Methane  is  an  attractive  fuel  gas  as  it  is  present  at  high  concentra¬ 
tions  in  several  fossil  and  renewable  sources  such  as  natural  gas 


[14,15]  and  biogas  [16].  The  utilization  of  methane  as  a  fuel  in  SOFC 
has  been  widely  studied  [17-20].  Methane  could  be  directly  used  as 
a  fuel  in  an  anode  chamber  of  SOFC  stack;  the  performance  of 
methane-fueled  SOFC  stack  was  significantly  lower  than  that  of  Re¬ 
fueled  SOFC  stack  as  described  in  our  recent  work  [21  ].  Moreover, 
direct  utilization  of  methane  in  SOFC  stack  caused  inhibition  of  the 
performance  of  SOFC  stack  due  to  carbon-deposited  blocking  inside 
anode  pore  [22].  Therefore,  the  pre-reformer  is  preferable  to  be 
installed  in  SOFC  system  in  order  to  initially  convert  methane  into 
hydrogen.  An  installation  of  pre-reformer  not  only  inhibits  carbon 
blockage  in  the  anode  chamber  but  also  enhances  the  power  density 
of  the  SOFC  stack  [23].  In  this  study,  the  average  stack  power  density 
is  considered  as  a  performance  indicator.  Operating  parameters  and 
SOFC  stack  configuration  are  varied  to  examine  operation  patterns, 
inlet  temperature  and  operating  voltage  that  offer  maximum  power 
density  under  stringent  thermal  constraints  i.e.  maximum  accept¬ 
able  stack  temperature  gradient  (Tc.acp)  and  maximum  acceptable 
stack  temperature  (Tmax,acp)-  Moreover,  a  new  configuration  of 
installing  inter-cooler  with  non-uniform  potential  SOFC  operation  is 
proposed  to  achieve  higher  average  power  density. 

2.  Modeling 

The  operation  of  SOFC  takes  place  at  the  following  configuration. 
The  reformed  methane  is  fed  into  the  SOFC  stack  at  the  entrance  of 
anode  chamber,  while  air,  performing  as  the  oxidizing  agent,  is  fed 
at  the  entrance  of  cathode  chamber  of  SOFC  stack.  Anode  and 
cathode  chambers  are  separated  from  each  other  by  solid  electro¬ 
lyte.  Inside  the  cathode,  oxygen  present  in  air  is  reduced  to  oxygen 
ions  (Eq.  (1))  which  permeate  via  the  solid  electrolyte  to  the  anode 
chamber  where  the  reaction  between  oxygen  ions  and  hydrogen 
takes  place  to  generate  steam  (Eq.  (2)).  CO  electro-oxidation  was 
assumed  to  be  negligible  due  to  its  lower  reaction  rate  compared 
with  H2  electro-oxidation  [24],  and  to  be  consumed  in  the  anode 
chamber  via  water-gas-shift  (WGS)  reaction,  where  reaction  rates 
are  high  at  elevated  operating  temperatures  [25-27].  The  SOFC  is 
composed  of  an  anode  chamber  (Ni-yttria  stabilized  zirconia),  an 
electrolyte  (yttria  stabilized  zirconia),  and  a  cathode  chamber 
(Lanthanum  strontium  manganite-yttria  stabilized  zirconia). 

lo2  +  2e-^02-  (1) 


h2  +  o2-^h2o 


(2) 


The  open  circuit  voltage  (0)  was  calculated  using  the  Nernst 
equation  (Eq.  (3)). 
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Fig.  1.  A  small  element  divided  for  the  calculation  in  SOFC  cell. 
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Irreversible  loss  generated  inside  the  SOFC  stack  was  calculated  in 
term  of  overpotential.  The  difference  between  open  circuit  voltage 
and  the  summation  of  all  irreversible  loss  was  named  operating 
voltage  (Eq.  (4)).  In  this  study,  three  major  types  of  overpotential 
were  taken  into  account:  activation  overpotential,  ohmic  over- 
potential  and  concentration  overpotential.  The  estimation  of  these 
overpotentials  was  given  in  our  recent  work  [28]. 


V  —  FJ  r/act  7/ohm  r/con 


(4) 


Inside  the  anode  chamber,  H2  was  consumed  via  the  electro¬ 
chemical  reaction.  Additional  hydrogen  is  obtained  via  methane 
steam  reforming  (Eq.  (5))  and  WGS  (Eq.  (6)).  The  rate  of  methane 
steam  reforming  reaction  was  computed  using  Eq.  (7)  [29],  while 
WGS  was  assumed  at  chemical  equilibrium  due  to  its  fast  reaction 
rates  at  high  operating  temperatures  [25-27].  Equations  of  mass  and 
energy  balances  were  solved  for  each  small  controlled  volume  (1- 
mm  distance)  (Fig.  1)  via  one-dimensional  (1-D)  finite  difference 
method  in  the  flow  direction  of  SOFC  stack  using  computer  coding  in 
Visual  Basic  for  MS  Excel.  The  following  assumptions  were  applied  in 
order  to  simplify  the  calculation:  (i)  SOFC  stack  was  operated  in 
adiabatic  mode;  (ii)  heat  radiation  between  solid  components  inside 
the  stack  was  neglected  [13];  (iii)  heat  conduction  in  the  solid 
electrolyte  was  neglected  [13];  (iv)  only  electrochemical  reaction  via 
H2  electro  oxidation  takes  place;  (v)  pressure  drop  generation  due  to 
the  flow  of  gas  along  the  stack  was  neglected.  The  calculation  of 
mass  balance  around  the  control  volume  was  performed  using  Eqs. 
(8)  and  (9)  for  the  anode  and  cathode  chambers,  respectively. 

CH4  +  H2O^CO  +  3H2  (5) 

C0  +  H20~H2+C02  (6) 
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Fig.  2.  Material  flow  schemes  of  a)  uniform  potential  SOFC  stack  (1 -section  SOFC),  and 
b)  non-uniform  potential  SOFC  with  two  sections  (2-section  SOFC). 


To  evaluate  the  viability  of  SOFC  stack,  two  parameters  were 
considered:  the  maximum  stack  temperature  gradient  ( TG )  and  the 
maximum  stack  temperature  (Tmax).  The  former  parameter  was 
defined  as  maximum  spatial  temperature  change  in  electrolyte 
along  the  stack  length  and  the  maximum  stack  temperature  was 
designated  as  maximum  value  of  temperature  along  the  SOFC  stack. 
The  maximum  acceptable  stack  temperature  gradient  (TG.acp)  of 
10  K  citT1  [30]  and  maximum  acceptable  stack  temperature 
(’fmax.acp)  of  1273  K  were  set  as  thermal  constraints  for  the  viable 
operation  of  SOFC  stack  in  this  study.  The  technical  terms,  i.e. 
average  current  density  (iavg),  average  power  density  (pav g),  total 
fuel  utilization  (l/f),  electrical  efficiency  (p)  and  %excess  air  (rair) 
were  defined  as  expressed  below; 


ref 


=  4274We 


(l  Pco(ph2)3  \ 

\  ^CH4^H20^ref/ 


<k  =<k- 


'  vk,nrn  ' 


vk,e\ec^^ 

2 T~ 
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vk,e\ec^Aj 


Km=Km+  2  F 


(9) 


Operating  voltage  was  assumed  to  be  constant  along  the  SOFC 
stack  as  its  current  collector  usually  has  a  high  electrical  conduc¬ 
tivity.  Heat  generation  in  electrolyte  due  to  the  irreversibility  of 
electrochemical  mechanism  and  exothermic  WGS  are  transferred 
into  the  bulk  gas  flow  inside  anode  and  cathode  chambers  through 
heat  convection  phenomena.  Energy  balance  around  solid  tri-layer 
for  each  control  volume  was  performed  using  Eq.  (10).  The  energy 
balance  of  bulk  fluid  flow  for  anode  and  cathode  chambers  was 
calculated  via  Eqs.  (11)  and  (12),  respectively. 


WT1  =  H>i  +  haAi(T{-T{) 


(11) 


Hf1  =  I-fc  +  hcAP  (t{  —  TJC) 
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Methane  was  assumed  to  be  initially  reformed  using  steam  to 
methane  molar  feed  ratio  of  2.5  operated  at  SOFC  feed  temperature 
(Tin)  prior  to  be  fed  to  the  anode  chamber  of  SOFC  stack  as  shown  in 
Fig.  2.  Chemical  equilibrium  was  assumed  for  this  reforming  reac¬ 
tion  in  the  pre-reformer.  The  term  “methane  flow  rate,  FCH4”  rep¬ 
resents  methane  feed  rate  to  the  reformer  and  the  term  “%  excess 
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Table  1 

Summary  of  model  parameters  used  in  1-D  analysis  [33—34]. 


Parameters 

Value 

SOFC  cell 

k 

750  pm 

lc 

50  pm 

L 

50  pm 

Cell  length 

400  mm 

Cell  width 

100  mm 

Anode  channel  height 

1  mm 

Cathode  channel  height 

1  mm 

ha  =  hc 

0.2  kj  m~2  s— 1  K1 

Operating  pressure 

1  bar 

air,  £:air”  was  defined  as  percent  of  oxygen  feed  rate  (to  cathode 
chamber)  to  the  stoichiometric  flow  rate  of  air  for  complete  com¬ 
bustion  of  Fch4-  The  values  of  fair  and  FCHa  were  kept  constant  at 
240%  and  1.08  mol  hr1,  respectively. 

3.  Results  and  discussion 

The  SOFC  stack  dimensions  and  input  parameters  were  initially 
set  as  summarized  in  Table  1.  Model  verification  for  the  calculation 
of  SOFC  stack’s  performance  was  performed  as  described  in  our 
previous  work  [31  ].  The  verifying  results  confirmed  that  the  uti¬ 
lized  SOFC  stack  model  could  well  predict  the  experimental  results 
of  Zhao  and  Virkar  [32]  for  the  mixture  of  97%  H2  and  3%  CO. 

To  improve  the  performance  and  operating  viability  of  SOFC,  the 
comparison  between  the  performance  of  1 -section  SOFC  (Fig.  2a) 
and  2-section  SOFC  (Fig.  2b)  was  first  investigated.  Inlet  tempera¬ 
ture  of  reformed  methane  and  air  was  set  at  998  K,  while  the 


operating  voltage  of  SOFC  stack  was  varied  until  its  optimum  value 
which  offers  the  highest  average  power  density  was  obtained.  As 
shown  in  Fig.  3a  and  b,  with  the  optimum  operating  voltage  of 
0.778  V,  maximum  power  density  was  found  at  the  stack  length  of 
0.25  m  while  stack  temperature  and  temperature  gradient  were 
always  lower  than  their  maximum  acceptable  values  along  the 
SOFC  stack.  Concentration  of  hydrogen  in  anode  gas  continuously 
decreased  in  the  flow  direction  (Fig.  3b)  due  to  the  utilization  of 
hydrogen  in  the  electrochemical  reaction.  In  the  first  half  of  SOFC 
stack,  the  power  density  increased  in  the  flow  direction  until  its 
maximum  value  was  reached  due  to  the  decrease  in  overpotential 
(increase  in  stack  temperature).  However,  near  the  gas  exit  of  the 
stack,  the  power  density  decreased  along  the  flow  direction.  This 
was  due  to  both  the  decrease  in  hydrogen  concentration  (Fig.  3b) 
and  rate  of  electrochemical  reaction.  Moreover,  the  open  circuit 
voltage  became  very  low  near  the  gas  exit  (Fig.  3c)  due  to  the  high 
stack  temperature.  This  would  imply  that  the  effect  of  the  decrease 
in  hydrogen  concentration  in  anode  bulk  gas  and  open  circuit 
voltage  dominated  the  effect  of  the  decrease  in  overpotential  near 
the  gas  exit  of  SOFC  stack.  From  the  calculations,  for  the  case  of  1  - 
section  SOFC,  the  optimum  operating  voltage  of  0.778  V  could 
provide  average  power  density  and  fuel  utilization  of  0.29  W  cm-2 
and  0.645,  respectively. 

From  the  results  observed  with  the  1 -section  SOFC,  the  differ¬ 
ence  between  the  open  circuit  voltage  and  operating  voltage  is  very 
high  and  becomes  lower  at  the  exit  point  of  the  stack  (Fig.  3c).  This 
envisions  different  behaviors  inside  the  SOFC  stack  and  suggests 
that  division  of  SOFC  stack  into  various  sections  with  different 
operating  voltage  may  improve  the  average  power  density  of  SOFC 
stack.  The  performance  of  2-section  SOFC  was  therefore  considered 
using  two  scenarios:  (a)  2-section  SOFC  with  high  voltage  at  the  1st 


Fig.  3.  The  distributions  of  a)  stack  temperature  and  power  density,  b)  stack  temperature  gradient  and  H2  concentration  in  anode  bulk  gas,  and  c)  open  circuit  voltage/operating 
voltage/overpotential  along  the  stack  for  1 -section  SOFC  at  maximum  average  power  density  ( Tin  =  998  K  with  optimal  operating  voltage  of  0.778  V). 
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section  (2S-1  V)  and  (b)  2-section  SOFC  with  high  voltage  at  the 
2nd  section  (2S-2  V).  For  the  section  operating  at  high  operating 
voltage,  the  operating  voltage  and  the  low-voltage  section  were  set 
to  0.8  V  and  0.76  V,  respectively.  To  evaluate  the  stack  performance 
for  these  scenarios,  the  1 -section  SOFC  operating  at  the  optimum 
operating  voltage  was  considered  as  a  base  case.  The  performance 
of  the  2-section  SOFC  illustrated  in  terms  of  power  density,  oper¬ 
ating  voltage  and  fuel  utilization  is  demonstrated  in  Figs.  4  and  5. 
For  2S-1  V  scenario  (Fig.  4),  the  power  density  at  the  1st  section  is 
lower  than  the  base  case  although  the  operating  voltage  was  higher 
than  that  of  the  base  case.  This  is  because  the  increase  in  temper¬ 
ature  and  also  fuel  utilization  were  lower,  and  the  power  density 
was  therefore  inhibited.  However,  2S-1  V  offered  a  higher  power 
density  at  the  2nd  stack  section  compared  with  the  base  case.  This 
could  be  explained  by  the  fact  that  there  were  large  amounts  of 
remaining  H2  fed  to  the  2nd  section  of  anode  chamber  due  to  low 
fuel  utilization  in  the  1st  section  of  stack,  leading  to  high  electro¬ 
chemical  reaction  rate  at  the  2nd  section  of  SOFC  stack  and  hence  a 
higher  power  density.  Due  to  high  fuel  utilization  at  the  2nd  section 
of  stack,  the  increase  in  stack  temperature  of  2S-1  V  was  also 
higher  when  compared  with  that  of  the  base  case.  For  2S— 2  V 
scenario  (Fig.  5),  superior  power  density  to  the  base  case  was 
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Fig.  5.  The  distribution  of  a)  power  density  and  stack  temperature  and  b)  open  circuit 
voitage/operating  voltage/overpotential  and  accumulated  fuel  utilization  along  the 
stack  between  the  base  case  and  2-section  SOFC  with  Vi  =  0.76  V  and  V2  =  0.8  V 
(7-in  =  998  I<). 


Stack  Length  /  in 


Fig.  4.  The  distribution  of  a)  power  density  and  stack  temperature  and  b)  open  circuit 
voltage/operating  voltage/overpotential  and  accumulated  fuel  utilization  along  the 
stack  between  the  base  case  and  2-section  SOFC  with  Vi  =  0.8  V  and  V2  =  0.76  V 
(Tin  =  998  K). 


Stack  Length  /  ni 

Fig.  6.  Comparison  of  the  distribution  of  stack  temperature  gradient  along  the  stack 
between  the  base  case,  2-section  SOFC  with  Vi  =  0.76  V  and  V2  =  0.8  V,  and  2-section 
SOFC  with  Vi  =  0.8  V  and  V2  -  0.76  V  (Tin  =  998  I<). 
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Fuel  Utilization 

Fig.  7.  Effect  of  operating  voltage  (Vi  and  V2)  on  the  average  power  density  and  fuel 
utilization  of  2-section  SOFC  (Tin  =  998  K). 


achieved  at  the  1  st  stack  section  due  to  higher  fuel  utilization.  This 
also  caused  a  sharp  increase  of  temperature  in  the  1st  stack  section 
for  2S-2  V.  However,  the  power  density  of  2S-2  V  in  the  2nd  stack 
section  became  lower  when  compared  with  that  of  the  base  case. 
This  might  be  due  to  lower  H2  concentration  and  rate  of  electro¬ 
chemical  reaction.  These  results  imply  that  the  difference  in  the 
open  circuit  voltage  and  operating  voltage  majorly  controlled  the 
fuel  utilization  as  shown  in  Fig.  5b,  the  voltage  difference  increased 
along  with  the  fuel  utilization  rate.  The  accumulated  fuel  utilization 
presented  in  Figs.  4b  and  5b  for  2S-1  V  and  2S-2  V,  respectively, 
indicates  the  total  fuel  utilization  from  the  gas  entry  of  SOFC  stack 
to  the  specified  stack  position.  An  increase  in  the  voltage  difference 
represents  a  high  overpotential  making  the  temperature  gradient 
at  the  1st  stack  section  for  2S-2  V  and  at  the  2nd  stack  section  for 
2S— 1  V  over  the  limit  as  shown  in  Fig.  6.  It  could  be  concluded 
that  the  decrease  in  operating  voltage  enhances  the  power  density; 
however,  thermal  constraint  of  SOFC  operation  could  also  be 
infringed.  By  considering  the  reason  above,  the  selection  of  sui¬ 
table  operating  voltage  for  2-section  stack  should  be  carefully 
considered. 

The  values  of  operating  voltage  for  the  1st  section  (Vi)  and 
2nd  section  (V2)  of  2-section  SOFC  stack  were  varied  to  study 
their  effects  on  average  power  density  and  fuel  utilization  as 
shown  in  Fig.  7.  The  results  demonstrate  that  the  decrease  in  V\ 
and  V2  could  improve  the  average  power  density  and  also  pro¬ 
mote  fuel  utilization  of  SOFC  stack.  However,  it  is  not  likely  that 
operation  of  2-section  SOFC  offered  a  greater  average  power 


b 


Stack  Length  /  ni 

Fig.  8.  The  distribution  of  a)  stack  temperature  and  power  density  and  b)  stack  tem¬ 
perature  gradient  along  the  stack  between  the  base  case  and  the  optimum  operation  of 
2-section  SOFC  (Tin  =  998  K). 

density  than  the  base  case  under  the  safe  operation  area  (the  area 
where  the  operating  conditions  offered  Tq  and  Tmax  below  their 
maximum  acceptable  values).  In  this  work,  the  optimum  values  of 
V\  and  V2  of  2-section  SOFC  were  found  to  be  at  0.768  and 
0.782  V,  respectively  (Fig.  7).  As  summarized  in  Table  2,  an 
insignificant  improvement  of  0.69%  in  average  power  density  of 
2-section  SOFC  over  1 -section  SOFC  is  observed.  At  the  optimum 


Table  2 

Comparison  of  the  performance  of  the  optimum  cases  of  1 -section  SOFC,  2-section  SOFC,  2-section  SOFC-T  and  3-section  SOFC-T. 


Operation  mode 

Optimum  case 
of  1 -section  SOFC 

Optimum  case  of 
2-section  SOFC 

Optimum  case  of 
2-section  SOFC-T 

Optimum  case  of 
3-section  SOFC-T 

FcH4,eq  (mol  S"1) 

0.0003 

0.0003 

0.0003 

0.0003 

Cair  (%) 

240 

240 

240 

240 

Uf(-) 

0.6451 

0.6462 

0.7739 

0.7425 

Operating  voltage  (V) 

0.778 

Vi  =  0.768 

Vi  =  0.79 

Vi  =  0.789 

V2  =  0.782 

V2  =  0.731 

V2  =  0.755 

V2  =  0.719 

Average  Power  density  (W  cm-2) 

0.2903 

0.2923 

0.3274 

0.3090 

Average  Current  density  (A  cm  2) 

0.3732 

0.3738 

0.4479 

0.4297 

ri  (%) 

39.13 

39.40 

44.13 

41.65 

Electricity  produced  (W) 

116.14 

116.93 

130.96 

123.60 

Maximum  stack  temperature  (I<) 

1271.71 

1273.11 

1227.57 

1169.14 

Maximum  stack  temperature  gradient  (K  cm-1) 

9.21 

9.93 

9.94 

9.97 
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Fig.  9.  Effect  of  operating  voltage  and  operating  temperature  on  power  density  of  SOFC 
calculated  based  on  the  chemical  equilibrium  composition  of  steam  methane 
reforming  reaction  with  the  steam  to  carbon  ratio  of  2.5. 

operation,  the  profiles  of  power  density,  stack  temperature,  and 
stack  temperature  gradient  along  the  SOFC  stack  of  2-section 
SOFC  were  almost  similar  to  those  of  the  base  case  as  shown  in 
Fig.  8a  and  b.  This  is  due  to  the  narrow  thermal-safe  operation 
area  which  limited  the  operation  of  2-section  SOFC  (Fig.  7).  It 
could  be  concluded  that  the  operation  of  2-section  SOFC  with  the 
consideration  of  a  solely  operating  voltage  cannot  significantly 
improve  the  SOFC  performance. 


a 


As  discussed  previously,  although  the  stack  temperature  played 
the  major  role  in  controlling  the  stack  performance,  thermal 
constraint  is  another  important  factor  that  could  limit  its  perfor¬ 
mance.  In  Fig.  9,  it  is  shown  the  effect  of  operating  voltage  and 
operating  temperature  on  the  power  density  of  SOFC  operation, 
which  were  calculated  based  on  the  chemical  equilibrium  compo¬ 
sition  of  methane  steam  reforming  reaction.  The  results  show  that 
the  power  density  was  improved  when  SOFC  was  operated  at  high 
temperature  and  low  operating  voltage.  To  further  envision  the 
effect  of  stack  temperature  on  the  SOFC  performance,  the  effect  of 
inlet  temperature  of  fuel  and  air  (Tm)  was  studied  for  1 -section 
SOFC.  The  results,  as  demonstrated  in  Fig.  10,  show  that  despite 
the  fact  that  the  operation  of  SOFC  at  high  temperature  and  low 
voltage  could  offer  a  high  power  density,  this  operation  was  still 
impractical  as  Tq  and  Tmax  exceeded  their  limiting  values.  As  shown 
in  Fig.  10a  and  b,  when  Tin  was  higher  than  998  K,  the  operation  of 
1 -section  SOFC  was  controlled  by  the  constraint  of  Tmax.acp,  while 
the  constraint  of  iG.acp  became  an  operation-controlling  criteria  as 
the  inlet  fuel  temperature  was  lower  than  998  K.  Moreover,  for  high 
Tin»  1  -section  SOFC  could  not  be  operated  at  low  operating  voltage 
due  to  the  thermal  constraints.  To  overcome  the  thermal  con¬ 
straints,  values  of  the  operating  voltage  larger  than  0.81  V  were 
required  when  Tin  is  1048  K,  while  a  value  of  the  operating  voltage 
as  low  as  0.778  V  could  be  operated  for  1 -section  SOFC  operating  at 
Tin  of  998  K  (Fig.  10a).  Operation  of  SOFC  at  a  higher  T[n  (ranging 
from  1048  I<  to  998  I<)  caused  lower  average  power  density 
(Fig.  10a)  and  total  fuel  utilization  (Fig.  10c). 

To  operate  at  favored  conditions  (high  temperature  and  low 
operating  voltage)  with  the  thermal-secured  operation,  the  new 
SOFC  configuration,  “non-uniform  potential  and  2-section  SOFC 
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Stack  Length  /  m 
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Fig.  10.  The  distribution  of  a)  stack  temperature  and  power  density,  b)  stack  temperature  gradient,  and  c)  accumulated  fuel  utilization  along  the  SOFC  stack  for  the  optimum 
operation  of  1 -section  SOFC  at  different  Tin. 
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Fig.  11.  Schematic  diagram  of  non-uniform  potential  2-section  SOFC  stack  equipped 
with  cooler,  2-section  SOFC-T. 
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Fig.  12.  The  distribution  of  stack  temperature  and  power  density  along  the  SOFC  stack 
for  the  optimum  case  of  2-section  SOFC-T  operating  at  different  Tin. 


stack  equipped  with  a  cooler,  (2-section  SOFC-T)”  was  proposed. 
The  scheme  of  this  configuration  was  similar  to  that  of  2-section 
SOFC,  except  installation  of  an  additional  cooler  between  the 
first  and  the  second  sections  of  the  stack  as  shown  in  Fig.  11.  With 
this  configuration,  the  outlet  gas  streams  from  the  cathode  and 
anode  chambers  of  the  1  st  section  of  stack  were  fed  to  the  cooler, 
where  temperature  of  gas  stream  was  reduced  to  7^,2  prior  to  be 
fed  to  the  2nd  section  of  stack.  The  effects  of  inlet  temperatures 
(Tinti)  and  (Tin,2)  of  fuel  and  air  for  the  1st  and  2nd  stack  sections, 
respectively,  as  well  as  those  of  operating  voltages  (Vi)  and  (V2)  of 
the  1  st  and  the  2nd  stack  sections,  respectively,  on  performance  of 
2-section  SOFC-T  were  studied.  To  simplify  the  calculations,  Tin>i 
was  assumed  equal  to  7^,2  and  both  of  them  were  similar  to  Tin. 


Fig.  13.  The  distribution  of  stack  temperature  gradient  along  the  SOFC  stack  for  the 
optimum  operation  of  2-section  SOFC-T  at  different  Tin. 


Fig.  14.  Power  density  and  optimum  V\t  and  V2  of  2-section  SOFC-T  at  different  T[n. 


The  operating  voltages  V\  and  V2  were  varied  until  their  optimum 
values  (the  value  that  allowed  the  maximum  average  power 
density  of  2-section  SOFC-T)  were  achieved  and  all  thermal  con¬ 
straints  were  pursued.  The  stack  temperature  and  power  density 
profiles  for  the  optimum  case  of  2-section  SOFC-T  along  the  stack 
with  the  variation  of  T[n  are  illustrated  in  Fig.  12.  The  results  show 
that  Tmax  of  2-section  SOFC-T  operated  at  Tin  between  998  K  and 
1098  K  were  largely  lower  than  those  of  the  base  case  (1 -section 
SOFC)  and  TmaX)acp.  Therefore,  the  thermal  constraint  for  the 
operation  of  2-section  SOFC-T  was  the  controlling  of  TG  of  SOFC 
stack  to  be  lower  than  TG>aCp.  As  shown  in  Fig.  13,  unlike  the  1- 
section  SOFC  and  2-section  SOFC,  the  limitation  of  the  value  of 
Tg  is  challenged  with  the  operation  of  2-section  SOFC-T,  especially 
for  the  operation  at  7in  =  1048  K,  where  this  situation  is  found  at 
both  gas  entry  and  middle  position  of  the  two  SOFC  sections.  The 
thermal  constraint  challenge  at  the  gas  entry  was  caused  by  the 
high  H2  concentration  of  anode  fuel  gas  stream  that  accelerates 
the  rate  of  electrochemical  reaction  at  this  stack  position,  while 
the  challenge  at  the  middle  of  the  stack  was  caused  by  the 
enhancement  in  the  electrochemical  reaction  rate  due  to  the  in¬ 
crease  in  stack  temperature.  Flowever,  the  selection  of  the 
appropriate  Tin  is  also  an  important  issue  for  the  control  of  the 
performance  of  2-section  SOFC-T.  At  7in  =  998  K,  the  thermal 
constraint  challenge  for  2-section  SOFC-T  is  found  only  at  near  the 
gas  exit  of  both  two  sections.  For  the  case  that  7in  =  1098  K,  the 
thermal  constraint  was  challenged  only  at  the  gas  entries  of  each 
section.  This  is  because  the  rate  of  electrochemical  reaction  at  the 
gas  entry  was  very  high  due  to  its  high  operating  temperature. 
Moreover,  with  high  7in,  desirable  low  operating  voltage  could  not 
be  operated  due  to  the  thermal  constraint  at  the  gas  entry.  The 
summary  of  average  power  density  and  optimum  operating 
voltage  at  different  Tin  of  2-section  SOFC-T  are  given  in  Fig.  14.  The 
maximum  power  density  of  0.327  W  cnrr2  was  achieved  at 
Tin  =  1048  K  and  optimum  values  of  V\  and  V2  of  0.79  and  0.731  V, 
respectively.  The  maximum  average  power  density  of  2-section 
SOFC-T  is  12.77%  higher  than  that  of  the  base  case  (1 -section 
SOFC).  If  Tin  was  lower  than  its  optimum  value,  average  power 
density  of  2-section  SOFC-T  became  lower  due  to  lower  electro¬ 
chemical  reaction  rate.  Increase  in  power  density  of  2-section 
SOFC-T  was  also  inhibited  when  operating  at  high  Tin  (Fig.  14). 
This  is  due  to  the  fact  that  low  desirable  operating  voltage  cannot 
be  operated  at  high  Tin  owing  to  the  thermal  constraint.  The  2- 
section  SOFC-T  could  offer  a  higher  stack  performance  compared 
with  the  1 -section  SOFC. 
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Fig.  15.  Schematic  diagram  of  non-uniform  potential  3-section  SOFC  stack  equipped  with  two  coolers,  3-section  SOFC-T. 


To  further  evaluate  the  performance  of  non-uniform  potential 
SOFC  stack  equipped  with  a  cooler,  the  non-uniform  potential  SOFC 
stack  with  3  sections  equipped  with  2  coolers,  3-section  SOFC-T, 
was  proposed  as  shown  in  Fig.  15.  With  this  configuration,  the 
length  of  each  stack  section  was  set  to  be  similar  to  one  another 
while  inlet  temperatures  of  fuel  gas  and  air  streams  for  each  stack 
section  were  assumed  similar  and  equal  to  T[n 
(T^  =  Tin,i  =  Tini2  =  Tin>3).  Tin  of  3-section  SOFC-T  was  varied  and 
optimum  values  of  V\,  V2 ,  and  V3  for  each  Tin  were  examined.  As 
shown  in  Fig.  16,  similar  to  the  case  of  2-section  SOFC-T,  the 
maximum  average  power  density  of  3-section  SOFC-T  was  achieved 
when  Tin  =  1048  K.  To  achieve  the  maximum  average  power  density 
for  each  Tm,  V\  was  always  set  to  be  higher  than  V2  and  V2  to  be 
higher  than  V3.  It  is  observed  that  the  increase  in  the  number  of 
sections  could  not  improve  the  SOFC  performance.  As  summarized 
in  Table  2,  maximum  achieved  average  power  density  of  3-section 
SOFC-T  was  0.309  W  cm-2,  which  was  lower  than  that  of  2-section 
SOFC-T  (0.327  W  cm-2).  As  shown  in  Fig.  17a,  since  the  length  of 
each  stack  section  of  3 -section  SOFC-T  was  too  short,  its  power 
density  could  not  reach  its  maximum  value.  With  the  defined  SOFC 
stack  dimension  and  electrochemical  reaction  rate,  unlike  2-section 
SOFC-T,  adequately  high  stack  temperature  could  not  be  achieved  at 
the  gas  exit  of  each  section  for  3-section  SOFC-T  as  shown  in 
Fig.  17b.  This  result  reveals  that,  apart  from  the  operating  voltage, 
average  power  density  as  well  strongly  depended  on  the  stack 
temperature. 

It  is  noted  that  this  work  presented  the  success  of  applying  non- 
uniform  potential  SOFC  with  a  novel  inter-cooler  to  overcome  the 
limitation  of  temperature  gradient  by  trial  methodology.  However, 
using  complete  optimization  methodology  could  offer  more  precise 
suitable  operating  parameters.  In  addition,  the  effects  of  other 
design  parameters  such  as  the  length  of  each  stack  section,  the  gas 
inlet  temperature  for  each  stack  section,  and  %  excess  air  should 
also  be  considered. 


Fig.  16.  Power  density  and  optimum  Vi,  V2 ,  and  V3  of  3-section  SOFC-T  at  different  Tin. 


4.  Conclusion 

In  this  study,  one-dimensional  finite  difference  method  was 
applied  to  evaluate  the  performance  of  non-uniform  potential  SOFC 
stack.  Power  density  was  considered  as  a  performance  indicator  of 
the  operation  of  SOFC  stack  fueled  by  reformed  methane.  The  eval¬ 
uation  was  performed  based  on  two  thermal  constraints:  maximum 
acceptable  stack  temperature  of  1273  K  and  maximum  acceptable 
stack  temperature  gradient  of  10  I<  cm-1.  Initially,  the  first  proposed 
non-uniform  potential  SOFC  with  two  sections  (2-section  SOFC)  was 
evaluated.  With  this  configuration,  SOFC  stack  was  divided  into  two 
sections  operating  at  different  operating  voltages.  The  values  of 


a 


b 


Fig.  17.  The  distribution  of  a)  power  density  and  b)  stack  temperature  along  stack 
length  of  the  optimum  cases  of  1 -section  SOFC,  2-section  SOFC-T  and  3-section  SOFC-T. 
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operating  voltage  of  each  section  of  2-section  SOFC  were  varied  to 
determine  maximum  average  power  density  under  the  two  thermal 
constraints.  The  results  indicated  that  2-section  SOFC  can  offer  only 
0.69%  average  power  density  higher  than  traditional  1  -section  SOFC 
with  the  optimum  values  of  operating  voltage  of  0.768  V  for  the  1  st 
section  and  of  0.782  V  for  the  2nd  section.  To  achieve  a  high  power 
density  and  simultaneously  overcome  the  thermal  constraints,  a  new 
configuration  of  non-uniform  potential  SOFC,  non-uniform  potential 
and  2-section  SOFC  stack  equipped  with  a  cooler  (2-section  SOFC-T) 
was  proposed.  In  this  configuration,  the  cooler  was  installed  to 
reduce  temperatures  of  the  fuel  gas  and  air  outlet  from  the  1  st  sec¬ 
tion  prior  to  be  fed  to  the  2nd  section.  It  was  found  that  not  only 
operating  voltage  that  played  an  important  role  in  the  operation  of  2- 
section  SOFC-T  but  also  the  inlet  fuel  and  air  temperatures  are  critical 
factors  affecting  the  performance  of  2 -section  SOFC-T.  The  value  of 
the  temperature  should,  therefore,  be  well  selected  to  overcome  the 
thermal  constraints  and  minimize  the  overpotential.  The  maximum 
average  power  density  of  the  2 -section  SOFC-T  was  achieved  when 
the  inlet  fuel  temperature  was  set  at  1048  K  and  the  values  of 
operating  voltages  of  the  1st  and  2nd  sections  were  equal  to  0.79  and 
0.731  V,  respectively.  By  operating  at  this  condition,  the  2-section 
SOFC-T  offered  12.8%  improvement  in  average  power  density 
compared  with  the  1 -section  SOFC.  With  the  defined  SOFC  stack 
dimension  and  electrochemical  reaction  rate,  an  increase  in  the 
number  of  section  from  2  sections  to  3  sections  (3 -section  SOFC-T) 
could  not  improve  the  average  power  density  since  the  high  stack 
temperature  was  not  reached  at  the  gas  exit  of  each  section  for  3- 
section  SOFC-T  and  its  overpotential  was  high.  Further  detailed 
study  should  be  performed  in  order  to  optimize  other  process  pa¬ 
rameters  of  2-section  SOFC-T,  such  as  the  length  of  each  stack  section 
and  the  fuel  gas  inlet  temperature  for  each  stack  section. 
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Nomenclature 

A:  SOFC  stack  active  area  [m2] 

E:  electromotive  force  [V] 

Eom.  electromotive  force  at  standard  pressure  [V] 

F:  Faraday  constant  (9.6495  x  104)  [C  mol-1] 

Fair, in:  inlet  air  flow  rate  to  SOFC  [mol  s_1] 

FcH4,eq:  methane  equivalent  flow  rate  [mol  s-1] 

Ffiiei'ini  inlet  fuel  flow  rate  to  SOFC  [mol  s-1] 
h:  heat  transfer  coefficient  [W  m-2  K-1] 

A heat  of  reaction  for  reaction  i  [J  mol-1] 

H:  energy  rate  [J  s-1] 
i:  current  density  [A  cm-2] 

I<ref.  equilibrium  constant  for  steam  reforming  reaction  [Pa2] 
l:  thickness  [pm] 

LHVfueiin.:  lower  heating  value  of  SOFC  inlet  fuel  [J  mol-1] 

Nj(:  molar  flow  rate  of  component  k  [mol  s-1] 
pi<:  partial  pressure  of  component  k  [Pa] 
p;  power  density  [W  cm-2] 

P:  electricity  produced  [W] 

rn:  reaction  rate  for  reaction  n  [mol  s-1] 

R:  gas  constant  (8.3145)  Q  mol-1  K-1] 

T:  temperature  [K] 

Tq:  maximum  stack  temperature  gradient  [K  cm-1] 

FG,acp:  maximum  acceptable  stack  temperature  gradient  [K  cm-1] 

Tin:  SOFC  inlet  temperature  of  fuel  and  air  streams  [K] 

Tmax'-  maximum  stack  temperature  [K] 

Tmax,acp:  maximum  acceptable  stack  temperature  [K] 

Uf, :  total  fuel  utilization  [— ] 

V:  operating  voltage  [V] 

Vn:  operating  voltage  of  section  n  [V] 

Greek  letters 

Eair-  %  excess  air  [%] 

r\:  electrical  efficiency  [%] 

Pact:  activation  loss  [V] 
ricon:  concentration  loss  [V] 

Pohm- :  ohmic  loss  [V] 

vicf.  stoichiometric  coefficient  of  component  k  in  reaction  i  [— ] 

Superscript 

j.  jth  controi  volume  (0,  1,  2,...,  400) 

Subscripts 
a:  anode  chamber 
avg:  average 
c:  cathode  chamber 
e:  electrolyte 

elec:  electrochemical  reaction 
k:  component  (CH4,  CO,  CO2,  H2,  H2O,  O2,  N2) 
n:  reaction  (ref,  WGS,  elec) 
ref:  reforming  reaction 
WGS:  water  gas  shift  reaction 


